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Abstract 

According to the Hohenberg-Kohn theorem, there is an invertible one-to-one relationship be- 
tween the Hamiltonian H oi a, system and the corresponding ground state density p(r). The 
extension of the theorem to the time-dependent case by Runge and Gross states that there is an 
invertible one-to-one relationship between the density p{rt) and the Hamiltonian H[t). In the proof 
of the theorem, Hamiltonians H/H{t) that differ by an additive constant C/ function C{t) are 
considered equivalent. Since the constant C / function C{t) is extrinsically additive, the physical 
system defined by these differing Hamiltonians H/H{t) is the same. Thus, according to the theo- 
rem, the density p(r) / p{rt) uniquely determines the physical system as defined by its Hamiltonian 
H/H{t). Hohenberg-Kohn, and by extension Runge and Gross, did not however consider the case 
of a set of degenerate Hamiltonians {H}/{H{t)} that differ by an intrinsic constant C/function 
C{t) but which represent different physical systems and yet possess the same density p{r)/p{rt) 
. The intrinsic constant C/function C(t) contains information about the different physical systems 
and helps differentiate between them. In such a case, the density p(r)/p{rt) cannot distinguish 
between these different Hamiltonians. In this paper we construct such a set of degenerate Hamil- 
tonians {H} / {H (t)} . Thus, although the proof of the Hohenberg-Kohn theorem is independent 
of whether the constant C/function C{t) is additive or intrinsic, the applicability of the theorem is 
restricted to excluding the case of the latter. The corollary is as follows: Degenerate Hamiltonian 
{H}/{H{t)} that represent different physical systems, but which differ by a constant C/function 
C{t) , and yet possess the same density p{r)/p{rt), cannot be distinguished on the basis of the 
Hohenberg-Kohn/Runge-Gross theorem. 
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I. INTRODUCTION AND COROLLARY 

This paper provides further insight into the first of the two Hohenberg-Kohn (HK) 
theorems [l| that constitute the rigorous mathematical basis for density functional 
theory. According to the theorem, for a system of electrons in an external field 
jrext^^^ _ -(-j-^g ground state electronic density p(r) for a nondegenerate state 

determines the external potential energy f (r) uniquely to within an unknown trivial 
additive constant C. Since the kinetic energy T and electronic-interaction potential 
energy JJ operators are known, the Hamiltonian H is explicitly known. The solutions \& 
of the corresponding time-independent Schrodinger equation, for both ground and excited 
states, then determine the properties of the electronic system. The wave function is thus a 
functional of the density: ^ = "^[p], and therefore all expectations are unique functionals of 
the density. Thus, the ground state density p(r) determines all the properties of the system. 

In the extension of the first HK theorem to the time-dependent case, Runge and 

n 

Gross(RG) [2] prove that for a system of electrons in a time- dependent external field 
T^^^{Yt) = — Vt'(rt), such that the potential energy v{rt) is Taylor- expandable about some 
initial time to, the density p{rt) evolving from some fixed initial state "^{to), determines the 
external potential energy uniquely to within an additive purely time- dependent func- 
tion C{t). Again, as the kinetic and electron- interaction potential energy operators are 
already defined, the Hamiltonian H{t) is known. The solution \E'(t) of the time- dependent 
Schrodinger equation then determines the system properties. Equivalently, the wave 
function \E'(t) is a functional of the density, unique to within a time- dependent phase. As 
such all expectation values are unique functionals of the density, the phase factor cancelling 
out. 

In the preamble to their proof, HK/RG consider Hamiltonians H/H(t) that differ by 
an additive constant C/function C(t) to be equivalent. In other words, the physical 
system under consideration as defined by the electronic Hamiltonian remains the same 
on addition of this constant/function which is arbitrary. Thus, measurement of properties 
of the system, other than for example the total energy E/E{t), remain invariant. The 
theorem then proves that each density p{r)/p{rt) is associated with one and only one 
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Hamiltonian H/H{t) or physical system: the density p{r)/p{rt) determines that unique 
Hamiltonian H/H{t) to within an additive constant C/function C{t). 

HK/RG, however, did not consider the case of a set of Hamiltonians {H}/{H{t)} that 
represent different physical systems which differ by an intrinsic constant C/function 
C{t), but which yet have the same density p{v)/p{vt). By intrinsic constant C/function 
C{t) we mean one that is inherent to the system and not extrinsically additive. Thus, this 
constant C/function C{t) helps distinguish between the different Hamiltonians in the set 
{H} / {H{t)}, and is consequently not arbitrary. That the physical systems are different 
could, of course, be confirmed by experiment. Further, the density p{r)/p{rt) would then 
not be able to distinguish between the different Hamiltonians {H}/{H{t)} or physical 
systems, as it is the same for all of them. 

In this paper we construct a set of model systems with different Hamiltonians 
{H}/{H{t)} that differ by a constant C/function C{t) but which all possess the same 
density p(r)/p(rt). This is the Hooke's species: atom, molecule, all positive molecular 
ions with number of nuclei J\f greater than two. The constants C/function C(t) contain 
information about the system, and are intrinsic to distinguishing between the different 
elements of the species. 

The corollary to the HK/RG theorem is as follows: Degenerate Hamiltonians {H} / {H{t)} 
that differ by a constant C/function C(t) but which represent different physical systems all 
possessing the same density p{v)/ p{vt) cannot be distinguished on the basis of the HK/RG 
theorem. That is, for such systems, the density p{v)/p{vt) cannot determine each external 
potential energy v{r)/v{rt), and hence each Hamiltonian of the set {H}/{H{t)}, uniquely. 

In the following sections, we describe the Hooke's species for the time-independent and 

time-dependent cases to prove the above corollary. 
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II. HOOKE'S SPECIES 

A. Time-Independent case. 

Prior to describing the Hooke's species, let us consider the following Coulomb species 
of two-electron systems and M nuclei: the Helium atom {M = 1; atomic number Z = 2), 
the Hydrogen molecule {M = 2; atomic number of each nuclei Z = 1), and the positive 
molecular ions {J\f > 2; atomic number of each nuclei Z — 1). 



In atomic units, the Hamiltonian of the Coulomb species is 

Hj^^f + U + V^ (1) 
where T is the kinetic energy operator: 

T--li:v^, (2) 
^ i=i 

U the electron-interaction potential energy operator: 



|ri - r2 

and V/v" the external potential energy operator: 



U = T-^, (3) 



VM = J2Mr^): (4) 



1=1 

with 



Mr) = E/c(r-R,). (5) 

where 

Here ri and r2 are positions of the electrons, — 1, ...jV) the positions of the nuclei, and 
fc{r — 'Rj) the Coulomb external potential energy function. Each element of the Coulomb 
species represents a different physical system. ( The species could be further generahzed 
by requiring each nuclei to have a different charge.) 
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Now suppose the ground state density p(r) of the Hydrogen molecule were known. Then, 
according to the HK theorem, this density uniquely determines the external potential energy 
operator to within an additive constant C: 

" ~|ri-Ri| ~ Iri-RsI ~ |r2-Ri| ~ R"^R^' 
Thus, the Hamiltonian of the Hydrogen molecule is exactly known from the ground state 
density. Note that in addition to the functional form of the external potential energy, the 
density also explicitly defines the positions Ri and R2 of the nuclei. 



The fact that the ground state density determines the external potential energy operator, 
and hence the Hamiltonian may be understood as follows. Integration of the density leads 
to the number of the electrons: / p{r)dr = N. The cusps in the electron density which 
satisfies the electron-nucleus cusp condition P], determine in turn the positions of the Af 
nuclei and their charge Z. Thus, the external potential energy operator Vj^ = Z]i^A/'(ri)) 
and therefore the Hamiltonian H are known. 



The Hooke's species comprise of two electrons coupled harmonically to a variable number 
Af of nuclei. The electrons are coupled to each nuclei with a different spring constants 
kj,j = 1...M. The species comprise of the Hooke's atom (A/" = 1, atomic number Z = 2, 
spring constant k), the Hooke's molecule {M = 2; atomic number of each nuclei Z = 1, spring 
constants ki and k2), and the Hooke's positive molecular ions (A/" > 2, atomic number of 
each nuclei Z = 1, spring constants ki,k2, ks...kj\/). The Hamiltonian Hj^ of this species is 
the same as that of the Coulomb species of Eq.(l) except that the external potential energy 
function is /ff(r — Rj), where 

/H(r-R,) = iA:,(r-R,)2. (8) 

Just as for the Coulomb species, each element of the Hooke's species represents a different 
physical system. Thus, for example, the Hamiltonian for Hooke's atom is 

Ha = - 1^2' + + \k[{v, - RO^ + {v, - (9) 

z z I Fx — r2 1 z 

and that of Hooke's molecule is 

gn^ = -^Vl^-^V2^+ , ^ , +i{fci[(ri-Ri)^ + (r2-Ri)Vfc2[(ri-R2)^ + (r2-R2)^]}, 
z Z \Y\ — r2 1 Z 

(10) 
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where k ^ k2, and so on for the various Hooke's positive molecular ions with J\f > 2. 



For the Hooke's species, however, the external potential energy operator Vj^ which is 



Vx^l E[fc,(ri - R,)^ + kj(v, - R,n (11) 



may be rewritten as 



Va/- = E - a)^ + (r2 - a)^] + C{{k}, {R},Af), (12) 



where the translation vector a is 



BL^j:k,R,/j:kj, (13) 

3=1 j=l 



and the constant C is 



with 



or 



C^b-d (14) 
b^J^k^Kl (15) 

d=(x:%R,)VE%, (16) 

3=1 3=1 

c = ^E'^.fc.(R.-R.)VEfcr (17) 

i,j=l j=l 

From Eq.(12) it is evident that the Hamiltonians H_\f of the Hooke's species are those of a 
Hooke's atom {Y.'jLikj — k), (to within a constant C{{k}, {R},jV)), whose center of mass 
is at a. The constant C which depends upon the spring constants {k}, the positions of the 
nuclei {R}, and the number J\f of the nuclei, differs from a trivial additive constant in that 
it is an intrinsic part of each Hamiltonian H^, and distinguishes between the different 
elements of the species. It does so because the constant C{{k}, {R},J\f) contains physical 
information about the system such as the positions {R} of the nuclei. 

Now according to the HK theorem, the ground state density determines the external 
potential energy, and hence the Hamiltonian, to within a constant. Since the density of 
each element of the Hooke's species is that of the Hooke's atom, it can only determine the 

6 



Hamiltonian of a Hooke's atom and not the constant C{{k}, {R}, jV). Therefore, it cannot 
determine the Hamiltonian Hj^ for jV > 1. This is reflected by the fact that the density of 
the elements of the Hooke's species does not satisfy the electron- nucleus cusp condition. ( 
It is emphasized that although the 'degenerate Hamiltonians' of the Hooke's species have 
a ground state wave function and density that corresponds to that of a Hooke's atom, 
each element of the species represents a different physical system. Thus, for example, a 
neutron diffraction experiment on the Hooke's molecule and Hooke's positive molecular ion 
would all give different results). 

It is also possible to construct a Hooke's species such that the density of each clement is 
the same. This is most readily seen for the case when the center of mass is moved to the 
origin of the coordinate system, i.e. for a = 0. This requires, from Eq.(13), the product of 
the spring constants and the coordinates of the nuclei satisfy the condition 

^%R, = 0, (18) 

i=i 

so that the external potential energy operator is then 

VAr) = lj:kjr' + lj:k,Kl (19) 

where r is the distance to the origin. If the sum Y^=\ kj is then adjusted to equal a particluar 
value of the spring constant k of Hooke's atom: 

E% = ^, (20) 

then the Hamiltonian Hj^ of any element of the species may be rewritten as 

Hj^{{k}, {R}, AT) = H,{k) + cm, {R}, AA), (21) 
where Ha{k) is the Hooke's atom Hamiltonian and the constant C{{k}, {R}, jV) is 

Ci{k},{R},Af)^j:kjRl (22) 

The solution of the Schrodinger equation and the corresponding density for each element 
of the species are therefore the same. 
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As an example, again consider the case of Hooke's atom and molecule. For Hooke's atom 
A/" = 1, Ri = 0, and let us assume k = ^. Thus, the external potential energy operator is 

Vair) = hr^ = V. (23) 

For this choice of k, the singlet ground state solution of the time-independent Schrodinger 
equation {Hj\f^ = -Eat^) is analytical 11^: 

^(rira) = De-^'/^e-^^'/^l + r/2), (24) 

where r = ri — r2,y = (ri + r2)/2, and D = 1/[2tt^^^{5^/tt + S)^^"^]- The corresponding 
ground state density p(r) = (^|p(r)|\E'), p(r) = '^('^ ^ ^i) is P, Ol 

p(r) = !^^D2g-rV2|7^ ^ ^3 ^ (8/V2^)re-"'/2 + 4(1 + r^)erf{r/V2}, (25) 

where 

er/(x) = ^ r e-'^dz. (26) 
V^r JO 

For the Hooke's molecule. A/" = 2,Ri = — R2, and we choose ki = k2 = |, so that the 
external potential energy operator is 

vUr) = lr' + j-^{Rl + Rl) = lr' + lR^ (27) 

where |Ri| = R. Thus, the Hamiltonian for Hooke's molecule differs from that of Hooke's 
atom by only the constant thereby leading to the same ground state wave function 
and density. However, the ground state energy of the two elements of the species differ by 

1 p2 

The above example demonstrating the equivalence of the density of the Hooke's atom 
and molecule is for a specific value of the spring constant k for which the wave function 
happens to be analytical. However, this conclusion is valid for arbitrary value of k for 
which solutions of the Schrodinger equation exist but are not necessarily analytical. For 
example, if we assume that for each element of the species (A/" > 2), all the spring constants 
kj,j = 1,2, ...A/" are the same and designated by k', then for the three values of k for the 
Hooke's atom corresponding to k = ^, |, 1, the values of k' for which the Hooke's molecule 
and molecular ion (A/" = 3) wave functions are the same are k' = k' — ^ ^- — ^ ^ 



8 



respectively. 

Thus, for the case where the elements of the Hooke's species are all made to have the 
same ground state density p(r), the density cannot, on the basis of the HK theorem, 
distinguish bewteen the different physical elements of the species. 

Corollary: Degenerate time-independent Hamiltonians {H} that represent different 
physical systems, but which differ by a constant C , and yet possess the same density p(r), 
cannot be distinguished on the basis of the Hohenberg-Kohn theorem. 
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B. Time-Dependent case. 



We next extend the above conclusions to the time-dependent HK theorem. Consider 
again the Hooke's species, but in this case let us assume that the positions of the nuclei are 
time-dependent, i.e. Rj = Rj(t). This could represent, for example, the zero point motion 
of the nuclei. For simplicity we consider the spring constant strength to be the same {k') 
for interaction with all the nuclei. The external potential energy vj^{rt) for an arbitrary 
member of the species which now is 



Mr^) = ^A:'5:(r-R,(t))^ (28) 



may then be rewritten as 



where at some initial time t^, we have Rj(to) = R-j.o- (Note that a spatially uniform time- 
dependent field F(t) interacting only with the electrons could be further incorporated by 
adding a term F(t) ■ r to the external potential energy expression.) The Hamiltonian of an 

element of the species governed by the number of nuclei M is then 

M 

Hj^{r,V2t) = Hm,o - k'J2[^jit) - %o] ■ (ri + r^) + C{k',Af,t), (30) 

i=i 

where Hj^fi is the time-independent Hooke's species Hamiltonian Eq.(21): 

H^fi = H^{k'), (31) 

and the time-dependent function 

C{k',Af,t) = k'j:[R'^{t)-Rl,]. (32) 

Note that the function C{k' ,M ,t) contains physical information about the system: in this 
case, about the motion of the nuclei about their equilibrium positions. It also differentiates 
between the different elements of the species. 

The solution of the time-dependent Schrodinger equation(if_A/(t)^'(t) = id'^{t)/dt) em- 
ploying the Harmonic Potential Theorem is 

^{T^Y^t) = exp{-t(l)it)}exp[~i{E^,ot - 2S{t) - 2^ • y}]^o(rr r^), (33) 
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where = - z{t),y = (ri + r2)/2, 

^^^^ = Ijr^^'^' " ^^^(^') V^', (34) 

the shift z{t) satisfies the classical harmonic oscillator equation 

M 

z{t) + kz{t) - k' ^[Ri(t) - Kj,o] = 0, (35) 

where the additional phase factor (f)(t) is due to the function C{k',M,t), 

(l){t) = [' C{k',Af,t')dt', (36) 

Jto 

and where at the initial time ^(rir2to) = which satisfies ^;v',o^o = -EV.o^o- Thus, 
the wave function ^(rir2t) is the time-independent solution shifted by a time-dependent 
function z{t), and multiplied by a phase factor. The explicit contribution of the function 
C{k',J\f,t) to this phase has been separated out. The phase factor cancels out in the 
determination of the density p{t) = {t)\f)\^ it)) = p(r — z{t)) which is the initial 
time-independent density p(rto) = Po(r) displaced by z{t). 

As in the time-independent case, the 'degenerate Hamiltonians' Hj^{viV2t) of the 
time-dependent Hooke's species can each be made to generate the same density p(ri) by 
adjusting the spring constant k' such that Hk' — k, and provided the density at the initial 
time to is the same. The latter is readily achieved as it constitues the time-independent 
Hooke's species case discussed previously. 

Thus, we have a set of Hamiltonians describing different physical systems but which can 
be made to generate the same density p{rt). These Hamiltonians differ by the function 
C{k',J\f, t) that contains information which differentiates between them. In such a case, the 
density p{rt) cannot distinguish between the different Hamiltonians. 

Corollary: Degenerate time-dependent Hamiltonians {H(t)} that represent different 
physical systems, but which differ by a purely time-dependent function C{t), and which 
all yield the same density p{rt), cannot be distinguished on the basis of the Runge-Gross 
theorem. 
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HOHENBERG-KOHN CASE 
■ff/ffW^/9(r)/p(rO 



|c/c(t) 1 

1 additive I 



v{t)/v{tt) - v'(T)/v'(it) = C/C{t) 



HOOKE'S SPECIES CASE 

.|c/C(t) |_{^}/{^(,)}^^(,)/^(,,) 
I intrinsic I 



FIG. 1: A schematic representation of the Hohenberg-Kohn theorem and its corollary. 
III. ENDNOTE 

The proof of the HK theorem is general in that it is vahd for arbitrary local form ( 
Goulombic^ Harmonic, Yukawa, oscillatory, etc.) of external potential energy v{r)/v{rt). (In 
the time-dependent case, there is the restriction that v{rt) must be Taylor-expandable about 
some initial time to.) For their proof, HK/RG considered the case of potential energies , and 
hence Hamiltonians, that differ by an additive constant C/function C{t) to be equivalent: 

v{r)/v{rt) - v'{r)/v'{rt) = C/C{t). (37) 

By equivalent is meant that the density p(r) / p{rt) is the same. The fact that the constant 
C/function C{t) is additive means that although the Hamiltonians differ, the physical 
system, however remains the same. The theorem then shows that there is a one-to-one 
correspondence between a physical system (as described by all these equivalent Hamil- 
tonians), and the corresponding density p{r)/p{rt). The relationship between the basic 
Hamiltonian H/H(t) describing a particular system and the density p{y)/ p{rt) is bijective 
or fully invertible. This case considered by HK/RG is shown schematically in Fig. 1 in 
which the invertibility is indicated by the double-headed arrow. 
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The case of a set of degenerate Hamiltonians {H}/{H{t)} that differ by a constant 
C / function C (t) that is intrinsic such that the Hamiltonians represent different physical 
systems while yet all possessing the same density p(r)/p(rt), was not considered by 
HK/RG. In such a case, the density cannot uniquely determine the Hamiltonian, and 
therefore cannot differentiate between the different physical systems. This case, also 
shown schematically in Fig.l, corresponds to the Hooke's species. The relationship between 
the set of Hamiltonians {H}/{H{t)} and the density p{r)/p{rt) which is not invertible is 
indicated by the single-headed arrow. 

We conclude by noting that the Hooke's species, in both the time-independent and time- 
dependent cases, does not contitute a counter example to the HK/RG theorem. The reason 
for this is that the proof of the HK/RG theorem is independent of whether the constant 
C/ function C{t) is additive or intrinsic. The Hamiltonians in either case still differ by a 
constant C/ function C{t). A counter example would be one in which Hamiltonians that 
differ by more than a constant C/ function C{t) have the same density p(r)/p(rt). 
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